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I. INTRODUCTION 

The concept of the strength of a system of partial differential equations (PDEs) was 
introduced by A. Einstein as a measure for the size of the solution space of such a system. 
In 6 A. Einstein defined the strength of a system of partial differential equations governing 
a physical field as follows: "... the system of equations is to be chosen so that the field 
quantities are determined as strongly as possible. In order to apply this principle, we propose 
a method which gives a measure of strength of an equation system. We expand the field 
variables, in the neighborhood of a point V, into a Taylor series (which presupposes the 
analytic character of the field); the coefficients of these series, which are the derivatives of 
the field variables at V, fall into sets according to the degree of differentiation. In every such 
degree there appear, for the first time, a set of coefficients which would be free for arbitrary 
choice if it were not that the field must satisfy a system of differential equations. Through 
this system of differential equations (and its derivatives with respect to the coordinates) the 
number of coefficients is restricted, so that in each degree a smaller number of coefficients 
is left free for arbitrary choice. The set of numbers of 'free' coefficients for all degrees of 
differentiation is then a measure of the 'weakness' of the system of equations, and through 
this, also of its 'strength'." 

Calculating by hand A. Einstein found out, that, for example, the potential and field 
formulations of Maxwell equations have different strengths for the dimension four. However, 
he did not obtain the exact expression of the above-mentioned number of free coefficients 
as a function of the degree of differentiation. Even though there were a number of works 
on the strength of a system of differential equations (in particular, on its relation to Cartan 
characters), see, for example, 20-22,25-28 , and 30 , there was no method of evaluating such a 
function until 1980 when A. Mikhalev and E. Pankratev 2,! showed that the strength of 
a system of algebraic partial differential equations (that is, a system of the form f\ = 
0, i G /, where are multivariate polynomials in unknown functions and their partial 
derivatives) is expressed by Kolchin's differential dimension polynomial associated with the 
differential field extension defined by the system. This observation allowed A. Mikhalev and 
E. Pankratev to develop two methods of determining the strength of a system of algebraic 
PDEs via computing the differential dimension polynomial of the corresponding differential 
field extension. The first method is based on construction of a characteristic set of the 
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ideal of differential polynomials defined by the system and then computing the differential 
dimension polynomial using the leading terms of the elements of the characteristic set (the 
idea of this approach comes from the original proof of Kolchin's theorem, see 9 (Chapter II, 
Theorem 6)). The second approach is based on the works by J. Johnson'' 8 , who showed that 
the differential dimension polynomial of a differential field extension can be computed as a 
Hilbert polynomial of the associated module of Kahler differentials. Using free resolutions 
for such a module, A. Mikhalev and E. Pankratev 23 evaluated the strength of several well- 
known systems of PDEs including the wave equation, both forms of Maxwell equations, 
Dirac equations (with zero mass), Lame equations, and some other systems of PDEs of 
mathematical physics. Note that A. Einstein, K. Mariwalla, M. Sue and some other authors 
who investigated the concept of strength in 1970s characterized the strength of a system by 
the "coefficient of freedom", an integer, that is fully determined by the leading coefficient 
of the differential dimension polynomial. The fact that such a polynomial provides a far 
more precise description of the strength than its leading term was justified by the result 
of W. Sit 2 " who proved that the set of differential dimension polynomials is well-ordered 
with respect to the natural order (/(£) < g(t) if and only if f(r) < g(r) for all sufficiently 
large integers r); this result allows one to distinguish two systems of PDEs with the same 
"coefficient of freedom" by their strength. 

Since 1980s the technique of dimension polynomials has been extended to the analy- 
sis of systems of algebraic difference and difference-differential equations. In a series of 
works whose results are summarized in 18 the second author proved the existence and devel- 
oped some methods of computation of dimension polynomials of difference field extensions 
and systems of algebraic difference equations. These polynomials determine A. Einstein's 
strength of a system of algebraic partial difference equations (we give the details in Section 
3 of this work) and, in particular, allow one to evaluate the quality of difference schemes for 
PDEs from the point of view of their strength. 

The next step in the analysis of systems of PDEs and systems of partial difference equa- 
tions is to consider their degrees of freedom with respect to different groups of basic operators 
(differentiations or translations). Theorems on multivariate dimension polynomials proved 
in 15 " 17 (see also 18 (Chapters 3, 4, 7)) allow to characterize the strength of a system of par- 
tial differential, difference or difference-differential equations in the case when the "weights" 
of basic operators of different groups are different. Methods of computations of multivari- 
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ate dimension polynomials for systems of differential, difference and difference-differential 
equations were developed i n 4 > 5 > 15 - 18 > 31 ; and 32 with the use of generalizations of the Grobner 
basis technique. In particular, the first author has implemented in Maple two algorithms of 
computation of bivariate difference-differential dimension polynomials via relative Grobner 
bases introduced in 31 . 

In this paper we present the theory and technique of differential, difference, and difference- 
differential dimension polynomials together with the applications of this technique to the 
analysis of fundamental systems of PDEs of mathematical physics and corresponding sys- 
tems of partial difference equations. In particular, we develop a method that allows one to 
compute the strength of such systems in the sense of A. Einstein and compare different dif- 
ference schemes for a given system of PDEs by their strength. We illustrate this method with 
the computation of the strength of the diffusion equation, Maxwell equations and equations 
for an electromagnetic field given by its potential, as well as with the computation of the 
strength of systems of difference equations obtained from these PDEs via different difference 
schemes. 

II. PRELIMINARIES 

In this section we present some basic concepts and results that are used throughout 
the paper. In what follows, N, Z, Q, and R denote the sets of all non-negative integers, 
integers, rational numbers, and real numbers, respectively. The number of elements of a set 
A is denoted by \A\. As usual, Q[t] denotes the ring of polynomials in one variable t with 
rational coefficients. By a ring we always mean an associative ring with unit element. Every 
ring homomorphism is unitary (maps unit element onto unit element), every subring of a 
ring contains the unit element of the ring. Unless otherwise indicated, by a module over a 
ring R we always mean a unitary left i?-module. 

2.1. Differential and difference rings and fields 

A differential ring (respectively, a difference ring) is a commutative ring R together with 
a finite set A = {S\, . . . ,S m } of mutually commuting mappings of R into itself such that 
each Si is a derivation of R (respectively, Si are injective endomorphisms of R also called 
translations) . The set A is said to be the basic set of the differential (or difference) ring R, 
which is also called a A-ring. If a A-ring is a field, it is called a A-field (this is a differential 
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field if A is a set of mutually commuting derivations and a difference field if the elements 
of A are endomorphisms) . If 6i,...,6 m are automorphisms of a difference ring R, we say 
that R is an inversive difference ring with the basic set A. In this case we denote the set 
{S±, ... , S m , 5± x , . . . , by A* and call R a A*-ring (if R is a field, it is called an inversive 
difference field or a A*-field). 

Let R be a A-ring and Rq a subring (ideal) of R such that 5(Rq) C Rq for any 5 G A. 
Then R is called a A-subring (respectively, a A-ideal) of R. If Rq is a A-subring of -R, we 
also say that R is a A-ring extension of .Ro- 
lf the elements of A act on R as mutually commuting derivations, we say that Rq is a 
differential subring (differential ideal) of R; if the elements of A are mutually commuting 
injective endomorphisms, we say that Rq is a difference subring (difference ideal) of R. 
Anyway, the prefix A-, depending on the context, means either "differential" or "difference", 
while the prefix A*- means "inversive difference". If R is a A*-ring (this assumption implies 
that A is a set of mutually commuting automorphisms of R), then a subring (ideal) Rq of R 
is called a A*-subring (respectively, A*-ideal) of R if a(Ro) C Rq for any a G A*. (A*-ideals 
are also called reflexive difference ideals of R; this term, as well as the term A*-ideal, is also 
used for ideals / of a difference A-ring R such that for any 5 G A, a G R, the inclusion 
5(a) G / implies a G I). If R is a A- (or A*-) field and Rq a subfield of R which is also a A- 
(respectively, A*-) subring of R, then R is said to be a A- (respectively, A*-) subfield of R 
while R is called a A- (respectively, A*-) field extension (or overfield) of Rq. In this case we 
also say that we have a A- (or A*-) field extension R/Rq. 

If R is a A-ring with a basic set A = {5±, . . . , 8 m }, then Ga (or if the set A is fixed) will 
denote the free commutative monoid generated by 6i, . . . , 8 m . Elements of G will be written 
in the multiplicative form 5^ . . . 5^ (k±, . . . , k m G N) and considered as the corresponding 
mappings of R into itself. If R is an inversive difference (A*-) ring, then Ta (or T if the set 
A is fixed) will denote the free commutative group generated by the set A. It is clear that 
elements of the group T (written in the multiplicative form 5^ . . . 5 l ™ where i±, . . . , i m G Z) 
act on R as automorphisms and is a subsemigroup of T. 

Let R be a A-ring and S C R. Then the intersection of all A-ideals of R containing 5* 
is denoted by [S]. Clearly, [S] is the smallest A-ideal of R containing S; as an ideal, it is 
generated by the set OS = {9(a) \6 G 9, a G S}. If J = [S], we say that the A-ideal J is 
generated by the set S called a set of A- generators of J. If S is finite, S = {a±, . . . , a^}, we 
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write J = [ai, . . . , a&] and say that J is a finitely generated A-ideal of the A-ring R. (In this 
case elements a±, . . . , are said to be IS.- generators of J.) If R is an inversive difference 
(A*-) ring and S C R, then the smallest A*-ideal of R containing S is denoted by [S]* (as an 
ideal, it is generated by the set TS = {7(a) |7 £ I\ a £ S 1 }. If S 1 is finite, S 1 = {ai, . . . , a*,}, 
we write [ai, . . . , a^]* for / = [S]* and say that / is a finitely generated A* -ideal of R; in this 
case the elements ax, . . . , a k are called A*-generators of I. 

Let R be a A-ring, i?o a A-subring of R and B O R. The intersection of all A-subrings 
of R containing R and B is called the A-subring of R generated by the set B over R ; it 
is denoted by Rq{B}. (As a ring, Rq{B} coincides with the ring Ro[{9(b)\b £ B, 9 £ 0}] 
obtained by adjoining the set {#(6)|6 £.£?,#£ 0} to the ring i? ). The set S is said to be 
the set of A-generators of the A-ring Rq{B} over Rq. If this set is finite, B = {b\, . . . , bk}, 
we say that R' = Rq{B} is a finitely generated A-ring extension (or A-ring extension) of 
Rq and write R' = Ro{bi, . . . If R is a A-field, Rq a A-subfield of R and i? C R, 

then the intersection of all A-subfields of R containing R and B is denoted by Rq{B) 
(or Ro(bi, . . . , bk) if B — {61, . . . , bk} is a finite set). This is the smallest A-subfield of R 
containing Rq and B\ it coincides with the field Ro({6(b)\b £ B,6 £ 0}). The set -B is 
called a set of A-generators of the A-field Rq(B) over i?o- If -R is an inversive difference 
(A*-) ring, Rq a A*-subring of R and B <Z R. Then the intersection of all A*-subrings of 
R containing Rq and B is the smallest A*-subring of R containing R and B. This ring 
coincides with the ring i? [{7(6)|6 £ B,j £ T}]; it is denoted by i?o{-E>}*- The set B is 
said to be a se£ 0/ A* -generators of i?o{-B}* over i? - If B — . . . , is a finite set, we 
say that S = Rq{B}* is a finitely generated A*-) ring extension (or A*-overring) of Rq and 
write S = Ro{h, . . . , bk}*. If -R is a A*-field, Rq a A*-subfield of R and B C R, then the 
intersection of all A*-subfields of R containing Rq and B is denoted by Rq(B)*. This field 
coincides with the field Ro({'y(b)\b £ -B, 7 £ T}). The set B is called a set of A* -generators of 
the A*-field extension R (B)* of Rq. If B is finite, = {&!,... , b k }, we write Ro{b±, . . . , b k )* 
for R {B}*. 

In what follows we often consider two or more A- (or A*-) rings R\, . . . ,R P with the 
same basic set A = {5\, . . . , 5 m }. Formally speaking, it means that for every i = 1, . . . ,p, 
there is some fixed mapping Vi from the set A into the set of all derivations or injective 
endomorphisms of the ring Ri such that any two mappings and Vi(5k) of Ri commute 

(1 < j, k < n). We shall identify elements 5j with their images Vi{8j) and say that elements 
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of the set A act as mutually commuting derivations or injective endomorphisms of the ring 
Ri (i = 1, . . . ,p). 

Let Ri and R2 be differential or difference or inversive difference rings with the same basic 
set A. A ring homomorphism : R\ — > R2 is called a A-homomorphism if 0(<5(a)) = <5(0(a)) 
for any 5 G A, a G R\. (Clearly, if : R\ — > R2 is a A-homomorphism of A*-rings, then 
<f)(5(a)) = 5(<f)(a)) for any 5 G A*, a G If R\ and i?2 are two A-overrings of the same 
A-ring R and : i?i — > R2 is a A-homomorphism such that 0(a) = a for any a G -Ro, we 
say that is a A-homomorphism over Rq or that leaves the ring R fixed. It is easy to see 
that the kernel of any A-homomorphism of A-rings : R — > R' is a A-ideal of R (moreover, 
in the case of difference rings, this kernel is a reflexive difference ideal of R\). Conversely, 
let g be a surjective homomorphism of a A-ring R onto a ring 5" such that Ker g is a A- or 
A*- (if R is a difference A-ring) ideal of R. Then there is a unique structure of a A-ring on 
5* such that g is a A-homomorphism. In particular, if J is a A- or A*- (if R is a difference 
A-ring) ideal of a A-ring R, then the factor ring R/I has a unique structure of a A-ring 
such that the canonical surjection R — > R/I is a A-homomorphism. In this case R/I is said 
to be the A-factor ring of R by /. 

If a A- (or A*-) ring R is an integral domain, then its quotient field Q(R) can be naturally 
considered as a A- (respectively, A*-) overring of R. (If A consists of derivations, then they 
extend to Q{R) via the quotient rule). In this case Q(R) is said to be the quotient A- 
(respectively, A*-) field of R. Clearly, if a A- (or A*-) field K contains an integral domain 
R as a A- (respectively, A*-) subring, then K contains the quotient A- (respectively, A*-) 
field Q(R). 

2.2. Differential, difference, and inversive difference polynomials. Algebraic 
differential and difference equations. 

With the above notation, let R be a A- (or A*-) ring and let U = {u\\\ G A} be a family 
of elements in some A- (respectively, A*-) ring extension of R. We say that the family U is 
A- algebraically dependent over R, if the family QU = {8(u\)\6 G 6, A G A} is algebraically 
dependent over R (that is, there exist elements v%, . . . ,Vk G QU and a non-zero polynomial 
f(Xi, . . . , Xk) with coefficients in R such that f(v\, . . . , Vk) = 0). Otherwise, the family U 
is said to be A- algebraically independent over R or a family of A-indeterminates over R. In 
the last case, the A-ring S = -R{(^a)aga} is called the algebra of A-polynomials over R. If A 
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consists of derivations (respectively, injective endomorphisms) , then S is also called a ring 
of differential (respectively, difference) polynomials in the difference (or a-) indeterminates 
{(^a)asa} over R. If R is a A*-ring and the family U considered above is A-algebraically 
independent over R, then the ring R{(u\)x^\}* is called the algebra of A "-polynomials in 
the A*-indeterminates u\ over R. 

If a family consisting of one element u is A-algebraically dependent over R, the element 
u is said to be A-algebraic over R. If the set {6(u)\6 E 0} is algebraically independent over 
R, we say that u is A-transcendental over the ring R. 

Let R be a A-field, L a A-field extension of R, and A E L. We say that the set A is 
A-algebraic over R if every element a £ A is A-algebraic over R. If every element of L is 
A-algebraic over R, we say that L is a A-algebraic field extension of i?. 

The following statement is proved in 9 (Chapter 1, Section 6), 3 (Chapter 2, Theorem I), 
and 10 (Propositions 3.3.7, 3.4.4) for differential, difference and inversive difference rings. 

Proposition II. 1 . Let R be a A- (respectively. A*-) ring and I an arbitrary set. Then 
there exists an algebra of A- (respectively, A*-) polynomials over R in a family of A- (respec- 
tively, A*-) indeterminates with indices from the set I. If S and S' are two such algebras, 
then there exists a A-isomorphism S —> S' that leaves the ring R fixed. If R is an integral 
domain, then any algebra of A- (respectively, A*-) polynomials over R is an integral domain. 

The algebra of A-polynomials over a A-ring R in a family of A-indeterminates with indices 
from a set J is a polynomial i?-algebra in the set of indeterminates QY = {y^jje/.ege with 
indices from the set / x 6. This algebra, as it is shown in 9 (Chapter 1, Section 6), 3 (Chapter 
2, Theorem I), and 111 (Propositions 3.3.7) can be viewed as a A-ring extension of R where 
HVi,d) = Vi,se for any 5 E A, y i>9 E BY. Setting y { = y iA we can write y Le as dy { . If R is a 
A*-ring, then the algebra of A*-polynomials over R in A-indeterminates with indices from 
a set J is a polynomial i?-algebra S in the set of indeterminates TY = {yi, 7 }ig/, 7 er with 
indices from the set I x T. As it is shown in (Propositions 3.4.4), S can be treated as a 
A* -ring extension of R where <5(?/j i7 ) — Ui,5^ 

for any 5 E A*, y irf E TY. In what follows we 
denote y iA by y { and write y ia (7 E T) as 7^. 

Let R be a A-ring, R{(yi)i e i} an algebra of difference polynomials in a family of A- 
indeterminates {(|/i)ig/}, and {(f]i) ieI } a set of elements in some A-ring extension of R. 
Since the set {6i)\i E 1,8 E 0} is algebraically independent over R, there exists a unique 
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ring homomorphism (fi v : R[(8yi) ie ifi e Q] — > -R[#(?7i)i6/,6»ee] that maps every 9 yi onto 9(rji) and 
leaves R fixed. Clearly, <p v is a surjective A-homomorphism of R{(yi) i€ j} onto R{(r]i) ie i}; 
it is called the substitution of (r/i)^/ for (yi)izi. Similarly, if R is a A*-ring, i?{(?/j)j 6 /}* an 
algebra of A*-polynomials over R and (?7j)jg/ a family of elements in a A*-ring extension of 
R, one can define a surjective A-homomorphism R{(yi)i e i}* — > R{(r]i)i e j}* that maps every 
yi onto r)i and leaves the ring R fixed. This homomorphism is also called the substitution of 
(Vi)iei f° r {Uijiei- (It will be always clear whether we talk about substitutions for difference 
(A-) or inversive difference (A*-) polynomials.) If g is a A- or A*- polynomial, then its image 
under a substitution of (r?i)i e / for (yi)iei is denoted by g((f]i) ie i). The kernel of a substitution 
is a A- (or A*- if we deal with difference or inversive difference polynomials) ideal of the 
A-ring R{(yi) i£ i} (respectively, of the A*-ring R{( yi ) i£ i}* if we consider substitution for 
inversive difference polynomials). This kernel is called the defining A- (or A*-) ideal of the 
family (r)i) ieI over R. 

If R is a A- (or A*-) field and {r]i)i & i is a family of elements in some A- (respectively, A*-) 
overfield S, then R{(r]i) ieI } (respectively, R{(i]i)i&i}*) is an integral domain (it is contained 
in the field S). It follows that the defining A- (or A*-) ideal P of the family (r^g/ over 
R is a prime A- (or A*- if we consider differences or inversive differences) ideal of the ring 
R{(yi)i(zi} (respectively, of the ring of A*-polynomials R{(yi)i £ i}*). Therefore, R{(r]i)i £ i) 
can be treated as the quotient A-field of the A-ring R{(yi)i & j}/P. (In the case of inversive 
difference rings, the A*-field R((r]i) ie j)* can be considered as a quotient A*-field of the 
A*-ring R{( yi ) ieI }*/P.) 

Let K be a A- (or A*-) field and s a positive integer. By an s-tuple over K we mean 
an s- dimensional vector a = (ai, . . . , a s ) whose coordinates belong to some A- (respectively, 
A*-) overfield of K. 

Definition II. 2 Let K be a A- (or A*-) field and let R be the algebra of A- (respectively, 
A*-) polynomials in finitely many A- (respectively. A*-) indeterminates yi, . . . ,y s over K . 
Furthermore, let $ = {fj\j G J} be a set of A- (respectively, A*-) polynomials in R. An s- 
tuple 7] — (771, . . . , rjs) over K is said to be a solution of the set $ or a solution of the system 
of algebraic A- (respectively, A*-) equations fj(yi, ■ ■ ■ ,y s ) — (j G J) if $ is contained in 
the kernel of the substitution of (rji, ... , r] s ) for (yi, ■ ■ ■ ,y s )- I n this case we also say that 77 
annuls $. 
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A system of algebraic difference equations $ is called prime if the A-ideal (or A*-ideal in 
the case of a system of difference or inversive difference equations) generated by <3> in the 
ring of A (or A*- if we deal with inversive difference equations) polynomials is prime. 

As we have seen, if one fixes an s-tuple rj = (771, . . . , rj s ) over a A- (or A*-) field K, then all 
A- (respectively, A*-) polynomials of the ring K{yi, . . . , y s } (respectively, K{y lt . . . , y s }*), 
for which rj is a solution, form a prime A- (respectively, A*-) ideal, the defining A- (respec- 
tively, A*-) ideal of 77. If $ is a subset of K{yi, ■ ■ ■ ,y s } (respectively, K{yi, . . . ,y s }*), then 
an s-tuple 77 = (771, . . . ,r) s ) over K is called a generic zero of $ if for any A- (respectively, 
A*-) polynomial /, the inclusion / G $ holds if and only if /(r/i, . . . , rj 8 ) = 0. 

Two s-tuples 77 = (771, . . . , rj s ) and £ = (Ci, • • • , Cs) over a A - (or A*-) field if are called 
equivalent over if if there is a A-homomorphism K(i]i, . . . , ?7 S ) — > if (G, . . . ,(s) (respectively, 
if (771, . . . , r] s )* — > if ((1, . . . , Cs)*) that maps each rji onto Q and leaves the field if fixed. 

Proposition II. 3 (see ( 3 (Chapter 2, Theorem VII), 10 (Propositions 3.2.6, 3.3.7)). Let R 
denote the algebra of A- (or A*-) polynomials in s A- (respectively, A*-) indeterminates 
. . . ,y s over a A- (respectively, A*-) field K . 

(i) A set $ ^ R has a generic zero if and only if $ is a prime A- ( or A*-, if we consider 
differences or inversive differences) ideal of R. If (771, . . . ,rj s ) is a generic zero of $, then 
K(r)i, . . . ,rj s ) is A-isomorphic to the A- (respectively, A*-) quotient field of R/<&. 

(ii) Any s-tuple over K is a generic zero of some prime A- (or A*-, if we deal with 
difference or inversive difference polynomials) ideal of R. If two s-tuples over K are generic 
zeros of the same prime A- (or A*-) ideal of R, then these s-tuples are equivalent. 

2.3. Ring of differential, difference, and inversive difference operators. Dif- 
ferential, difference, and inversive difference modules. 

Let R be a differential or difference ring with a basic set A = {5\, . . . , 5 m } and let O be 
the free commutative semigroup generated by A. If 6 = 5* 1 . . . 5^ e 6 (fci, . . . , k m G N), 
then the number ordO = Yl™=i is called the order of 9. Furthermore, for any r e N, the 
set {6 G G I ordO < r} is denoted by 0(r). 

Definition II. 4 An expression of the form Yleee ae @> where ae G R for any 9 G and only 
finitely many elements ae are different from ; is called a A-operator over the ring R. (If A 
is the set of mutually commuting derivations, then a A-operator is also called a differential 
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operator; if A consists of mutually commuting injective endomorphisms, a A-operator is 
called a difference operator ). Two A-operators J2eee a $ an< ^ Seee^ are considered to be 
equal if and only if ag = bg for all 9 G 0. 

The set of all A-operators over a A-ring R can be equipped with a ring structure if we 

J2eee ae ^^ 1 )' 5a = a5 + 5(a) (respectively, 5a = 5 (a) 6 if R is a difference ring and A is 
the basic set of endomorphisms of R) for any A-operators X^ee ag9, X^ee bg9 and for any 
a G R, 5 G A, and extend the multiplication by distributivity. The ring obtained in this way 
is called the ring of A-operators over R; it will be denoted by T>. (If A is a set of derivations, 
T> is also said to be the ring of differential operators over the differential ring R; if A is a 
set of endomorphisms, T> is called the ring of difference operators over R.) 

The order of a nonzero A-operator A = Yleee ae ® ^ ^ is defined as the number ordA = 
ma,x{ord9 \ ag ^ 0}. We also set ordO = — oo. 

Let T> r = {A G T>\ordA < r} for any r G N and let T> r = for any r G Z, r < 0. Then 
the ring T> can be treated as a filtered ring with the ascending filtration (T> r ) re z- Below, 
while considering T> as a filtered ring, we always mean this filtration. 

Definition II. 5 Let R be a A-ring and T> the ring of A-operators over R. Then a left T>- 
module is called a A-R-module. (If A is a set of derivations, we also use the term differential 
R-module; if A is a set of endomorphisms, we use the term difference R-module). In other 
words, an R-module M is a A-R-module if the elements of A act on M in such a way that 
5(x + y) = 5(x) + 5(y), 5(5'x) = 5\5x), and 5(ax) = 5(a)x + a5(x) (if A consists of deriva- 
tions, so A- means "differential") or 5(ax) = 5(a)5(x) (if A consists of endomorphisms, 
so A- means "difference") for any x, y G M; 5, 5' G A; a G R. If R is a A-field, then a 
A-R-module M is also called a vector A-R-space. 

If R is an inversive difference ring with basic set A = {5i, . . . ,5 m } and T is the free 
commutative group generated by A, then the order of an element 7 = S^ 1 . . . 5^ G T 
(ki, . . . , k m G Z) is defined as ore? 7 = Y^=i IM- Also, for any r G N, we set T(r) = {7 G 
T I ord^i < r}. In this case, by a A*-operator we mean an expression of the form Yl-y^r a il-> 
where a 7 G R for any 7 G V and only finitely many elements a 7 are different from 0. As in 
the case of A-operators, two A*-operators are considered to be equal if and only if for any 
7 G T, the coefficients of 7 in these operators are the same. 
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Clearly, a A*-ring R can be also treated as a A-ring, so every A-operator over R can be 
also considered as a A*-operator. The set of all A*-operators over the ring R can be naturally 
considered as a ring extension of the ring T> of A-operators over R where the operation are 
defined in the same way as they are defined in T> with additional rules 5~ 1 a = 5 _1 (a)5 _1 
and (W -1 = <5 _1 5 = 1 (a G R, 5 G A) extended by distributivity. The resulting ring will be 
denoted by V*; it is called the ring of A*-operators (or inversive difference operators) over 
R. A left P*-module is called a A*-i?-module (or an inversive difference -R-module). Such a 
module is actually a A- .R-module M with an additional action of the elements of the form 
d^ 1 (8 G A) such that 5(5~ 1 (x)) = 5~ 1 (5(x)) for any x G M (the other rules are the same 
as in Definition II. 5 except for that the elements S and 5' are taken from the set A* rather 
than from A). 

The order of a nonzero A*-operator A = Xl 7 er ^ ^* * s defined as the number 
oral A = ma,x{ord'~f | a 7 ^ 0}, and we also set ordO = — oo. The ring T>* will be treated as a 
filtered ring with the ascending filtration (T>*) r£ z where D* = {A G T>*\ordA < r} for any 
r G N and V* = for any r G Z, r < 0. 

III. DIFFERENTIAL AND DIFFERENCE DIMENSION POLYNOMIALS 

In this section we present main theorems on dimension polynomials of differential and 
difference modules and field extensions. Then we show how one can determine the strength 
of a system of partial differential or difference equations by computing the corresponding 
dimension polynomial. 

With the above notation, let R be a A- (respectively, A*-) ring. We say that a A-R- 
module (respectively, A*-i?-module) M is finitely generated, if it is finitely generated as a 
left T>- (respectively, T>*-) module. By a filtered A- (respectively, A*-) module we always 
mean a left T>- (respectively, T>*-) module M equipped with an exhaustive and separated 
filtration. Thus, a filtration of M is an ascending chain (M r ) r& i of -R-submodules of M such 
that V r M s C M r+S (respectively, V*M S C M r+S ) for all r, s G Z, M r = for all sufficiently 
small r G Z, and Urez M r = M. A filtration (M r ) r£ z of M is called excellent if all R- 
modules M r (r G Z) are finitely generated and there exists ro G Z such that M r = D r _ ro M ro 
(respectively, M r = V*_ rQ M ro ) for any r G Z, r > r . Note that if R is a A-field and M is a 
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finitely generated A- R- module, 



s 



i=i 



for some fi, ■ ■ ■ , f s £ M, then 



( 



i=i 



) 



is an excellent filtration of M, and a similar remark can be made about a finitely generated 
A*-i?-module. 

The following result combines theorems on dimension polynomials of differential and 
difference modules obtained in' and 12 (see also 10 (Theorems 5.1.11, 6.2.5 and Propositions 



Theorem III.l Let R be a A-field whose basic set consists of m operators (derivations or 
injective endomorphisms). Let V be the ring of A- operators over R, and let (M r ) r& i be an 
excellent filtration of a A-R-module M. Then there exists a polynomial ip(t) G Q[t] with the 
following properties. 

(i) ij}{f) = diHifj(M r ) for all sufficiently large r G Z, that is, there exists r G Z such that 
the last equality holds for all integers r > r . (as usual, dim^(M r ) denotes the dimension of 
the vector i?-space M r ). 

(ii) degtp(t) < m and the polynomial ip(t) can be written as 



(this polynomial takes integer values for all sufficiently large integer values of t). 

(iii) The integers d = degip(t), c m and Cd (if d < m) do not depend on the choice of the 
excellent filtration of M. Furthermore, c m is equal to the maximal number of elements of M 
linearly independent over the ring T>. 

The polynomial ip(t) whose existence is established by Theorem III.l is called the A- 
dimension polynomial (differential or difference dimension polynomial depending on the 
nature of the set A) of the A-i?-module M associated with the excellent filtration (M r ) r& z. 



5.2.12, 6.2.17)). 




rn 



where Cq, c±, . . . , c m G Z and 
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The integ ca are called the A-type, A- dimension, and typical A- dimension of 

M, respectively. A number of results on differential and difference dimension polynomials, 
as well as some methods of their computation, can be found in 1 " (Chapters 5-9). 

The following is an analog of Theorem III. 1 for inversive difference modules (see 10 (The- 
orem 6.3.3 and Proposition 6.3.15) or 18 (Theorems 3.5.2, 3.5.8)). 

Theorem III. 2 Let R be a A* -field whose basic set A consists of m automorphisms of 
R, and let (M r ) re z be an excellent filtration of a A* -R-module M . Then there exists a 
polynomial x(t) in one variable t with rational coefficients such that 

(i) x( r ) = dim^(M r ) for all sufficiently large r G Z; 

(ii) degx(t) < tti and the polynomial x(t) can be represented in the form 



where a G Z and o(t n ) is a polynomial in Q[t] of degree less than m. 

(iii) The integers a, d = degx{t) and the coefficient of t d in the polynomial x(t) do not 
depend on the choice of the excellent filtration of M. Furthermore, a is equal to the maximal 
number of elements of M linearly independent over the ring T>* . 

The polynomial x(t) is called the A* -dimension polynomial of the A*-i?-module M asso- 
ciated with the excellent filtration (M r ) re g- 

The next result combines Kolchin's theorem on differential dimension polynomial" (Chap- 
ter II, Theorem 6) and the corresponding result for difference field extensions proved in 12 . 

Theorem III. 3 Let K be a A-field whose basic set consists of m operators (derivations or 
endomorphisms) . Let L = K(r]i, . . . ,rj s ) be a A-field extension of K generated by a finite 
family 7] = {rji, . . . ,r] s }. Then there exists a polynomial 4>n\K{t) G Q[t] with the following 
properties. 

(i) <t>n\K{. r ) = trdegKK({9r]j \ 9 G 0(r), 1 < j < s}) for all sufficiently large r G N. 

(ii) deg4>r)\K{t) < n an d the polynomial (fi^it) can be written as 



x(t) 



a 



t rn + o{t m ) 



ml 




where a , . . . , a m G Z. 
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(hi) The integers a m ,d = deg (j} v \K(t) and are invariants of the polynomial (p^xit) , that 
is, they do not depend on the choice of a system of a -generators T). Furthermore, a m = A- 
trdegxL where A-trdegKL denotes the A-transcendence degree of L over K , that is, the 
maximal number of elements £1, • ■ • G L such that the family {9^i\9 G 6, 1 < % < k} is 
algebraically independent over K . 

The polynomial 4> v \K{t) whose existence is established by Theorem III. 3 is called the A- 
(differential or difference depending on the nature of A) dimension polynomial of the A-field 
extension L of K associated with the system of A- generators rj. The integers d = deg <fr v \K(t) 
and ad are called, respectively, the A-type and typical A-transcendence degree of L over K. 
These invariants of ^ v \k{P) are denoted by A-type K L and A-t.trdegxL, respectively. 

Notice that if the elements r]i, . . . ,rj s are A-algebraically independent over K (that is, 
the set {6rji \ 6 G 6, 1 < i < s} is algebraically independent over K) and <f> v \K{t) is the 
corresponding difference dimension polynomial of L/K (we use the notation of the last 
theorem), then 



for all sufficiently large r G N (CardO(r) is the number of solutions (ki, . . . , k m ) G W n of 
the inequality k\ + • • • + k m < r; it is well known (see, for example, 1 " (Proposition 2.1.9)) 



The following theorem shows the existence of a dimension polynomial of a finitely gener- 
ated inversive difference field extension. 

Theorem III. 4 Let K be a A*-field whose basic set A consists of m automorphisms of 
K . As before, let V be the free commutative group generated by A, and for any r G N, let 
T(r) = {7 G r I ord'-f < r}. Furthermore, let L = K(rji, . . . , t] s )* be a A* -field extension of K 
generated by a finite family rj = {rji, . . . , ?7 S }. Then there exists a polynomial ip n \K{t) G Q[t] 
with the following properties. 

(i) ip n \ K {r) = trdegKK({'yrjj \ 7 G T(r), 1 < j < s}) for all sufficiently large r G N. 



0r?|x( r ) = trdegKK({rrjj \ r G T, 1 < j < s}) = s ■ CardQ(r) = s 




that this number is ( r+m )). Therefore, in this case 
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(ii) deg ip v \K{t) < m and the polynomial ^^(t) can be written as 

^\K(t) = —t m + 0(t m ) 

ml 

where a G Z and o(t m ) is a polynomial in Q[t] of degree less than m. 

(hi) The integers a, d = deg i/j v \K(t) and the coefficient oft d in the polynomial iprj\K{t) do 
not depend on the choice of a system of generators rj. Furthermore, a = A-trdegxL. 

(iv) If 7)i, ... , 7] s are A- algebraically independent over K , then 

Let us consider a prime system of algebraic A- (differential or difference) or A*-(inversive 
difference) equations 

Ai( yi ,...,y s ) = Q (i = l,...,p) (1) 

where Ai(yi, . . . ,y s ) are A- (or A*- ) polynomials in the ring R = K{yi, . . . , y s } (respectively, 
in R = K{yi, . . . , y s }*) and let P be a prime A-ideal (respectively a prime A*-ideal if we 
consider the difference or inversive difference case) of R generated by the right-hand sides 
of system (1). Furthermore, let 7]i be the canonical image of ?/j in the factor ring R/P 
(1 < i < s). It is easy to see that for every r G N the intersection P(lR r is a prime ideal of the 
ring R r and the quotient helds of the rings R r /P PI R r and K[{9(r]j) \ 9 G 6(r), 1 < j < s}] 
(respectively, K\\p/(r)j) |7 G T(r), 1 < j < s}]) are isomorphic. Considering the case of 
algebraic differential or difference equations we can apply Theorem III. 3 and obtain that 
there exists a polynomial (pp(t) in one variable t with rational coefficients such that 

cj) P (t) = trdegicKftdfa) | 9 G 0(r), 1 < j < s}) = trdeg K (R r /P n Rr) 

for all sufficiently large r G Z, deg<pp(t) < m and the polynomial <f>p(t) can be written as 

m 
i=0 

where a , . . . , a m G Z and a m = A-trdegxK \{9(i)j) \ 9 G 9(r), 1 < j < s}). 

In the case of a system of difference equations (including the case when such a system 
involves negative degrees of basic translations, which act as automorphisms), one can apply 
Theorem III. 4 that shows the existence of a polynomial ipp{t) G Q[t] such that 

^ P (r) = trdegKKd^r/j) | 7 G T(r), 1 < j < s}) = trdeg K (R r /P n i? r ) 
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for all sufficiently large r £ Z, degip(t) < m and the polynomial ipp{t) can be written as 

^p(t) = -^* B + o(* B ) 

where ap = a-trdegxiR/P)- 

With the above notation, the numerical polynomial <pp(t) (respectively, ipp{t)) is called 
a differential (respectively, difference) dimension polynomial of system (1). It is also said to 
be a A- (respectively. A*-) dimension polynomial of the system. 

Taking into account Einstein's approach described in the Introduction, one can say that 
for all sufficiently large r, the value <pp(r) of the differential dimension polynomial <f>p(t) of a 
system of algebraic differential equations is the number of Taylor coefficients of order < r of 
an analytic solution that can be chosen arbitrarily. (These Taylor coefficients are the values 
of derivatives of order < r of the solution computed at the point in whose neighborhood 
we consider its expansion. The dependence of coefficients is understood as their algebraic 
dependence over the field of coefficients of the system.) Thus, 4>p(t) can be viewed as a 
measure of strength of system (1), so the problem of computation of differential dimension 
polynomials is important not only for the study of differential algebraic structures, but also 
for the study of equations of mathematical physics. 

Considering a system of equations in finite differences over a field of functions in several 
real variables, one can use Einstein's approach to define the concept of strength of such a 
system as follows (cf. Einstein's description of the strength of a system of PDEs presented 
in the Introduction). Let 

•U7: /.) (i = l,...,p) (2) 

be a system of equations in finite differences with respect to s unknown grid functions 
fx, . . . , f s in n real variables xx, . ■ . ,x n with coefficients in some functional field K. We also 
assume that the difference grid, whose nodes form the domain of considered functions, has 
equal cells of dimension hx x • • • X h n (hx, ■ ■ ■ , h n £ R) and fills the whole space IR n . As an 
example, one can consider a field K consisting of a zero function and fractions of the form 
u/v where u and v are grid functions defined almost everywhere and vanishing at at most 
finitely many nodes. (As usual, we say that a grid function is defined almost everywhere if 
there are at most finitely many nodes where it is not defined.) 
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Let us fix some node V and say that a node Q has order i (with respect to V) if the 
shortest path from V to Q along the edges of the grid consists of i steps (by a step we mean 
a path from a node of the grid to a neighboring node along the edge between these two 
nodes). Say, the orders of the nodes in the two-dimensional follows (a number 

near a node shows the order of this node). Let us consider the values of the unknown grid 
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FIG. 1: 2-dimensional grid 

functions /i, . . . , f s at the nodes whose order does not exceed r (r G N). If /i, . . . , f s should 
not satisfy any system of equations (or any other condition), their values at nodes of any 
order can be chosen arbitrarily. Because of the system in finite differences (and equations 
obtained from the equations of the system by transformations of the form fj{x±, . . . , x s ) t— > 
fj(xi + kihi, . . . , x s + k n h n ) with ki, . . . , k n G Z, 1 < j < s), the number of independent 
values of the functions fi, ■ ■ ■ , f s at the nodes of order < r decreases. This number, which 
is a function of r, is considered as a "measure of strength" of the system in finite differences 
(in the sense of Einstein). We denote it by S r . 

With the above conventions, suppose that the transformations <x,- of the field of coefficients 
K defined by 

(1 < j < n) are automorphisms of this field. Then K can be considered as an inversive dif- 
ference field with the basic set a = {a±, . . . , a n }. The replacement of the unknown functions 
fi by difference indeterminates (z = l,...,s) leads to a system of algebraic difference 
equations of the form (1). If this system is prime (e.g., we deal with a system of linear 
difference equations), then its difference dimension polynomial ip{t) expresses the strength 
S r . Thus, this polynomial can be naturally viewed as the measure of Einstein's strength of a 
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given system of equations in finite differences. In what follows, the A*-dimension polynomial 
ijj(t) will be called the difference dimension polynomial of the system. 

Methods of computation of A- and A*- dimension polynomials of a system of algebraic 
partial differential or difference equations developed so far are based either on building of a 
characteristic set of the considered above associated A- (or A*- ) ideal P in K{yi, . . . ,y s } 
or on constructing a free resolution of the module of Kahler differentials associated with the 
extension K(rji, . . . ,rj s ) (or K(r]i, . . . , rj s )*). The corresponding computations can be found, 
for example, in 10 (Chapter 9) and 18 (Chapter 7). The main drawback of the mentioned 
approaches is the lack of efficient algorithms for constructing characteristic sets and serious 
restrictions on the systems to which one can apply the method of free resolutions. In the 
last case, a system of difference equations with inversive difference operators is supposed 
to be linear and symmetric, that is, whenever an equation involves a A*-operator u = 

ai^ 11 . . . 8^ m H h a r 5i rl . . . 8^ m (a* G K), which contains a term a8^ ...S 1 ^ (a E K, a ^ 

0), then it also contains all terms of the form bdf l1 . . . 8 m lm with nonzero coefficients b E K 
and all 2 m distinct combinations of signs before li,...,l m . In what follows we explain a 
method of computation of dimension polynomials (and therefore, the strength of a system of 
algebraic partial differential or difference equations), which does not have these restrictions. 

Implementations for computing Grobner bases in modules of differential and difference 
operators are available, e.g., in the Mgfun package" for Maple or in the Plural extension of 
Singular 11 . 

IV. COMPUTATION OF THE STRENGTH OF A SYSTEM OF 
DIFFERENCE EQUATIONS VIA GROBNER AND GENERALIZED 
GROBNER BASIS TECHNIQUES. EXAMPLES 

Let K be a difference or inversive difference field of characteristic with basic set A = 
{5\, . . . , 8m}. As we have seen, the ring of A-operators over K carries many properties of a 
polynomial ring in m variables over K. In order to underline the relationship between the 
Grobner basis method for A-.ff-modules considered below and the classical Grobner basis 
technique for polynomial ideals we will denote the ring of A-operators over K by K[A] and 
set 

[A] :={8^---8t\k 1: ... : k m EN}. 
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Similarly, if A is a family of mutually commuting automorphisms of K, we set 

[A*] :={6?---6%> | h,...,k m eZ} 

and denote the ring of A*-operators over K by if [A*]. Furthermore, a free left if [A]- 
(respectively, if [A*]-) module with a set of free generators E = {ei, . . . , e q } will be denoted 
by K[A]E (respectively, K[A*]E) and [A]E (respectively, [A*]E) will denote the set of all 
elements of the form Ae^ where 1 < % < q and A G [A] (respectively, A G [A*]. Such elements 
of the free module are called terms. 

Let E = {ei, . . . , e q } be a finite set of free generators of the left if [A*] module K[A*]E 
and F C K[A*]E finite. There are two popular approaches for computing a Grobner basis 
of the left if [A*]-module k[a*]{F). The first approach is due to the second author 16 ' 18 
the idea also appears 19 - works by introducing new variables for the inverses &f , . . . , 5" 1 of 
81, . . . ,5 m and doing computations in the resulting free module of difference operators. The 
second approach, originated by the second author 13 , was enhanced by Winkler and Zhou 31 ' 32 
who introduced the concept of so-called generalized term orders therefore making K[A*]E a 
well-ordered set. In the following we will outline the first approach. Proofs for termination 
and correctness of the algorithms can be found, e.g, in 16,18 ' 19 . 

A. Computing Grobner bases of inversive difference modules via standard 
bases of associated difference modules 

Let if be a differential (respectively, difference) field with basic set A = {8%, . . . , d~ m } of 
derivations (respectively, endomorphisms) of if. 

Definition IV. 1 Let -< be a total order on the set of terms [A]E such that for all elements 
1 A, T), n G [A] , e, e' G E we have 

1. e -< Ae, 

2. fiXe -< firje' whenever Ae -< rje' . 

Then the relation -< is called an admissible order. 

For any / = ai/i + . . . + a n f n G K[A]E with a 1} . . . , a n G K, f h . . . , f n G [A}E and for a 
given admissible order -< we denote the highest term appearing in / with nonzero coefficient 
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Algorithm 1: Reduction_algorithm 

IN: / / e K[A]E, finite G C K[A]E, and an admissible order -<, 

OUT: r € K[A]E such that / is reducible to r modulo G and r is irreducible modulo G. 
r := f 

while there exist g € G and A £ [A] such that lt(Ag) = lt(r) do 
r -=r - Aj^G 
end while 
return r 

by lt(/), i.e., 

lt(/) := max{/j | 1 < % < n, a, 7^ 0}, 

and call it the leading term of /. The corresponding coefficient is called the leading coefficient 
of / and is denoted by lc(/). 

Definition IV. 2 Let f,g £ A'[A]£'\{0} and -< an admissible order. If there exists A £ [A] 
Wrt't/i lt(Ag) = lt(/) we say that f is reducible to h := / — modulo g and write 

Otherwise we say that f is irreducible modulo g. Let G = {g±, . . . ,g p } C iT[A]S \ {0}. If 
there exist n £ M, / , . . . , f n , i\, . . . , i n £ {1, . . . ,p} such that 

f = fo ^ffjj fl ^g i2 " ' ' ^g in fn—'-f 

we say that f is reducible to r modulo G. Otherwise we say that f is irreducible modulo G. 

The process of reduction is described by the following algorithm. 

Definition IV. 3 Let -< be an admissible order, N a submodule of K\E\E and G C N\ {0} 
finite such that every ^ / £ N is reducible to modulo G. Then G is called a Grobner 
basis of the module N. 

Every finitely generated A'[A]-module M has a Grobner basis that can be computed, e.g., 
via Buchberger's algorithm starting with any finite generating set G of M (see Algorithm 2 
below). 
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Algorithm 2: Buchberger ' s_algorithm 

IN: G C K[A]E \ {0} finite, -< an admissible order, 
OUT: G C if[A]£ \ {0} being a Grobner basis of K[A] {G). 
G:=G 

while there exists g,g' € G such that S(g,g') is not reducible to modulo G do 
G := G U {Reduction_algorithm(S'(5, 5'), G, -<;)} 
end while 
return G 

Definition IV. 4 Let -< be an admissible order on [A}E, gi,gi G K[A}E \ {0}, Ai,Ai G 
[A],ei,e2 G -E snc/i £/ia£ lt^(gi) = Aiei and lt_< (^72) = A 2 e2- The least common multiple 
lcm(lt x (£/i),lt^(p 2 )) o/lt^(^i) and lt^(# 2 ) «s defined by 



lcm(lt^(^),lt x (^)) := 



lcm(Ai, A 2 )ei if e\ = e 2 , 
if e 1 7^ e 2 . 



Let Mi,n 2 G [A] &e given 6?/ 



_ lcrnjlt^yt^ga)) _ lcm(lt^( qi ), lt^(g 2 )) 

Ul and M2 : " MS) 

Taen the S-polynomial S{gi,g2) of g± and g 2 is defined by 

o [91,92) ■= ui- — — — - - n 2 ; 



lc^(gx) lc^(g 2 ) 

The following theorem being a special case of 15 (Thm. 4.12) and 18 (Thm. 3.3.15) describes 
how to obtain the dimension polynomial associated with a system of differential or difference 
equations. 

Theorem IV. 5 Let M be a difference K-vector space generated (as a left K[A]-module) 
by elements mi, . . . , m q , F a free K[A]-module with set of free generators E = {ei, . . . , e q }, 
7r : F — > M the difference epimorphism (e, H- for i = 1, . . . , q) and N := ker(7r). Let 
G C K[A}E be a Grobner basis of N with respect to the term order -< defined by 

°1 °rn e i ^ °1 °m ^3 

* ' ' * {kl + ■ ■ ■ + k m , i, ki, . . . , km) <lcx (^1 + 1 ■ ' + lyni jiWi ■ ■ ■ i ^m); 
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where <\ ex denotes the lexicographic order. For r G Z let 
M r := {Am G [A] {mi, . . . , m q } | ord A < r} and 

U r := {Ae G [A]E | ord A < r and Ae ^ ^^(g) for any g G G, /i G [A]}. 

T/ien (M r ) rg 2 an excellent filtration of M and for any r G N £/ie set ir(U r ) is a basis for 
the K-vector space M r 

The following proposition is obtained from 1 " (Prop. 2.2.11.) by realizing that a term Ae is 
irreducible if and only if there exist no r] G [A], g G G with lt(g) = fie and rjp = A. 

Proposition IV. 6 With the notation of Theorem IV. 5 for every i = 1, . . . , q, let 

Gi := {lt( ) | g G G, \t(g) G [A]e,} 
and let A, = (Xij.k) £ pjl G il x2m satisfy the following condition: for every 

lt(s) = «r--<^---/t"e l GG l 

i/iere exzste j G {1, . . . , wit/i (A^x, . . . , Aj J)2m ) = (a 1; . . . , a m , &i, . . . , b m ). 

Furthermore, for any l,n G N with 1 < n and < I < n, let A(l,n) denote the set of 
all l-element subsets of {1, ... ,n} and for every 1 < i < q, 0^£G A(l, \Gi\), let \i£ t k '■ = 
max je5 Aij ifc and X it ^ k := 0. Finally, for any 1 < % < g,£ G iei fa := J2k^i\(,k- 

Then 

p^eb-d' e ( r+ tr /M )- 

An idea of the following kind was also by Ziming Li and Min Wu 19 . Let K be an inversive 
difference field with basic set of automorphisms A = {Si, . . . , S m } and F C K[A*]E a finite 
set of generators for a left if[A*] module. Let a±, . . . , a m , (3i, . . . , (3 m be endomorphisms of 
K such that for i = 1 , . . . , m and k G K we have 

ci!i(/c) := 5j(/c), and := ^(k). 

Then A' can be considered as a difference field with basic set E := {ai, . . . , a m , . . . , 
(3 m }- By p : AT[A*] — >■ A"[S]/x[s]({cKiA e — e I 1 < ? < m,e G A 1 }) we denote the natural 
isomorphism 

p : 5* 1 . . . 5^e ^ a^* 1 '"} . . . ^{^^^^{^b.o} ^{-fc^o^ (e £ ^ 
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Let F := p(F) U {a^e - e | 1 < i < m, e G E}. Then, K[A*]E / k[a*](F) is isomorphic 
to K[T]E / k[t,](F) , so in order to compute the A*-dimension polynomial of a A*-i^-module 
K[A*]E / k[a*]{F) associated with a finite system of generators F, it suffices to compute 
the E-dimension polynomial of the E-i^module K[H]E / k\t,]{F) associated with the set of 
generators F. 

B. Examples for the computation of differential and difference dimension 
polynomials 

In this subsection we give several examples for the computation of differential dimension 
polynomials associated with systems of differential equations arising in mathematical physics 
and of difference dimension polynomials associated with their difference schemes. 

Example IV. 7 (Diffusion equation in 1-space) 

The diffusion equation in one spatial dimension for a constant collective diffusion coefficient 
a and unknown function u(x,t) describing the density of the diffusing material at given 
position x and time t is given by 

du(x,t) d 2 u(x,t) 

~^t~ = a dx* ■ (3) 

Differential dimension polynomial 

Let K be a differential field with basic set A = {S x = J^,5 t = J)-} containing a and let 
M be a differential K -vector space generated as K[A]-module by one generator m satisfying 
the defining equation 

5 t m = ab 2 x m. 

Then M is isomorphic to the factor module of a free K[A]-module with free generator e by 
its submodule N generated by 

G := {5 t e — ad^e}. 

Since G consists of only one element there are no S-polynomials. Therefore G is already a 
Grobner basis of N for any admissible order on [A]e. Let the admissible order -< on [A]e be 
given by 

S^5^e -< S^S^e :<^> (k x + k t , k x , k t ) < lex (/ 
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Then for all 2 < r G N we have 

U r = {5 x x 5^e I k x + h<r, 5 x x 5^e is irreducible modulo G} 

= {e, 5 t e, 8 r t e, 5 x e, 5 x 5 t e, SJ^e}, 

and therefore \U r \ — 2r + 1. 

Thus, the differential dimension polynomial associated with the diffusion equation in one 
spatial dimension for a constant collective diffusion coefficient is given by <fi(r) — 2r + 1. 
Difference dimension polynomial for forward difference scheme 

In order to obtain a forward difference scheme for the diffusion equation (3) every oc- 
curence of 9U q x '^ and is replaced by u(x + l,t) — u(x,t) and u(x,t + 1) — u(x,t), 
respectively. We obtain 

u(x, t + 1) - u(x, t) = a(u(x + 2, t) - 2u(x + 1, t) + u(x, t)). (4) 

Let K be an inversive difference field with basic set A = {5 X : x > x + l,6t ■ t > t + 1} 
containing a and let M be an inversive difference K -vector space generated as a left K[A*]- 
module by one generator m satisfying the defining equation 

btvn — m = a(5 x m — 25 x m + m). 

Then M is isomorphic to the factor module of a free K[A*]-module with free generator e by 
its submodule N generated by 

G := {5 t e — a5^e + 2a5 x e — (1 + a)e}. 

We will compute the difference dimension polynomial associated with the difference scheme 
(4) using the method described at the end of subsection A. Thus, we consider K as a difference 
field with basic set E = {a x : x i— > x + 1, at : 1 1— > t + 1, fi x : x i— > x — 1, flt ". t h- > t — 1}. Let 
G := {gi := a t e — aale + 2aa x e — (1 +a)e, #2 := ct x f3 x e — e, gs := a t (3 t e — e} and I =a"[s] (G). 
Then K[T,]e/I is isomorphic K[A*]e/N via the isomorphism 

We fix an admissible order -< on [E]e defined by 
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(a x + a t + b x + b t , a x , a t , b x , b t ) < lex (c x + c t + d x + d t , c x , c t , d x , d t ) 

and compute a Grobner basis of I with respect to -<. The S-polynomial S{g\,g<i) of gi and 
g 2 is given by 

S{gi,g2) = -2a x /3 x e - -a t (3 x e+ ( 1 + - J x e+a x e 

a \ a J 

and is reducible modulo g 2 to — \at(5 x e + (l + -) f3 x e + a x e — 2e, which is irreducible modulo 
G. Hence, g± := — ^a t (3 x e+ (l + -j (3 x e + a x e — 2e should be inserted into G. For S(gi, #3) = 
-2a x a t f3 t e - \a 2 t f3 t e + (l + \) a t /3 t e + a 2 x e we have 

S '(91,93) — > 93 ~-a 2 f3 t e + a 2 x e+ ( 1 + - J a t /3 t e - 2a x e 



Furthermore 



>„ 3 ale + ( 1 + - J a t /3 t e - 2a x e - -a t e 
a J a 



r 9i 



> 93 0. 



1 + - J a t f3 t e- ( 1 + - }e 
a 1 \ a 



S(gi, 94) = aale + (a + l)alj3 x e — 2a x a t j3 x e a 2 j3 x e — 2aa 2 e 

a 

+ (l + ~J a t (3 x e 
— > 91 (a + l)a 2 x f3 x e - 2a x a t (5 x e - -a 2 t (5 x e + a x a t e 

(X 

+ f 1 + ~J a tPx£ - (1 + a)a x e 

— > gi -2a x a t p x e - -a 2 (5 x e + a x a t e + (2a + 2)a x (3 x e 
a 

+ (2 + 2/a)a t fi x e - (1 + aja^e + (-2 - a - l/a)(3 x e 



1 



r 92 



a t /3 x e + a x a t e + (2a + 2)a x /3 x e + (2 + 2/a)a t /3 x e 



a 

— (1 + a)a x e — 2a t e + (—2 — a — l/a)(3 x e 
gi (2a + 2)a 3t f} x e + f 1 + a^e - (1 + a)a x e 
+ (-2 - a - l/a)/3 x e 



r 92 



1 + - J a^e - (1 + a)a x .e + (-2 - a - l/a)/3 x e 
aj 

+ (2a + 2)e 



— >*0 

S{g2,9z) = a x (3 x e - a t (3 t e 
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r 92 





a t /3 t e + e 



r 93 

S(g2, 9a) = ac? x e + (a + l)a x (3 x e — 2aa x e — a t e 
— > gi (a + l)a x (3 x e - (1 + a)e 

S (#3, P4) = aa x fi t e + (a + l)/3 x /3 t e - (3 x e - 2a/3 t e 
=■ 95 

S(gi,gs) = - [1 + -) a x (3 x f3 t e + -a x f3 x e a t (5 t e 

\ a J a a 

+ (l + i) fte 

— ^92 -«x/3xe - -a t (3 t e 
a a 

— > q2 — a t /3 t e + -e 
a a 

5(^2,^) = -(l + -)(3 2 x f3 t e+-f3 2 x e + 2(3 x f3 t e-/3 t e 
\ a J a 

=■■ 96 

Further computations of the S -polynomials S(gi, g§), 1 < % < 5, show that all of them are 
reducible to modulo {gi, . . . ,g^}. Hence, a Grobner basis for I is given by 

gi = a t e — aa 2 x e + 2aa x e — (1 + a)e, 

g2 = a x p x e - e, 
#3 = a t/ 9 t e - e, 

#4 = —a t f3 x e + + f3 x e + a x e - 2e, 
a \ a J 

g 5 = aa x (3 t e + (a + l)/3 x f} t e - (3 x e - 2af3 t e, 
ge = ~ ( 1 + - J + -/3 x 2 e + 2/^e - fte 



a / a 
Applying Proposition IV. 6 we obtain 

| {A G [S]e | ord(A) < r, A is irreducible modulo G}\ = 5r. 

for all sufficiently large r. Therefore, the inversive difference dimension polynomial associ- 
ated with the difference scheme (4) is given by <f>{r) = 5r. 
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Difference dimension polynomial for symmetric difference scheme 

In order to obtain a space symmetric difference scheme for the diffusion equation (3) 

d 2 u(x,t) . du(x,t) , , , . N , . . . 

every occurrence of — — - — and — is replaced by u{x + 1, t) — 2u[x, t) + u{x — 1, t) 

and u(x, t + 1) — u(x, t), respectively. We obtain 

u(x, t + 1) — u(x, t) = a(u(x + 1, t) — 2u(x, t) + u(x — 1, t)). (5) 

Let K be an inversive difference field with basic set A = {5 X : x i— >■ x + 1, 8f : t t + 1} 
containing a and let M be an inversive difference K -vector space generated as a left K[A*]- 
module by one generator m satisfying the defining equation 

5 t m — m = a(5 x m — 2m + <5^ 1 m). 

Then M is isomorphic to the factor module of a free K[A*] -module with free generator e by 
its submodule N generated by 

G := {5 t e - a5 x e - 5~ l e + (2a - l)e}. 

Now consider K as a difference field with basic set £ = {a x > x + l,ctt : 1 1- > t + 1, (3 X : 
x (->• x-l,/3 t : 1 1-)- Let G := {g 1 := a t e - a x e - (3 x e + (2a - l)e, g 2 := a x (3 x e-e,g 3 := 
OitfitZ — e} and I =k[z] (G) . Then K[T]e/I is isomorphic K[A*]e/N via the isomorphism 

Let us fix an admissible order -< on [E]e defined by 

(a x + a t + b x + b u a X: a t , b x , b t ) < [ex (c x + c t + d x + d t , c x , c t , d x , d t ) 
A Grobner basis of I is then given by 

{gi := af3 2 x f3 t e - (2a - l)(3 x (3 t e - (3 X + a/3 t e, 

g 2 ■= --a t (3 x e + (3 x e- [ 2 - - ) (5 x e + e, 
a \ a J 

g 3 := a t (3 t e - e, 

:= aa x e — a t e + a(3 x e — (2a — l)e}. 

Applying Proposition IV. 6 we obtain 

| {A G [T]e | ord(A) < r, A is irreducible modulo G}\ = 4r 
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for all sufficiently large r. Hence, the inversive difference dimension polynomial associated 
with the difference scheme (5) is given by <p(r) = 4r. 

Thus, the symmetric difference scheme for the diffusion equation has higher strength 
(that is, smaller dimension polynomial) than the forward scheme, so the symmetric scheme 
is more preferable from the point of view of strength. 

Example IV. 8 (Maxwell equations for vanishing free current density and free 
charge density) 

Let E = (Ei, E 2 , E 3 ), D = (D h D 2 ,D 3 ),H = (Hi,H 2 ,H 3 ), B = (B h B 2 ,B 3 ), J f = 
(Ji, J21 J3) and pf be functions of (x,y,z,t) that denote electric field strength, electric 
displacement vector, magnetic field strength, magnetic displacement vector, free current 
density and free charge density, respectively. With 

v ._ ( 9 9 9 

\dx' dy' dz 

Maxwell's equations in 3 spatial dimensions are given by 

OB dD 
V-D = p f , V-B = 0, Vx£ + - = 0, and VxH=J f + —. 

at ot 

Assuming Jf = and pj = 0, Maxwell's equations can be considered as a set of homogeneous 
linear differential equations. 

Differential dimension polynomial 

Let K be a differential field with basic set A = {5 X — -^,6 y — -§^,S Z — -§^,8t — Jj}- 
Assuming Jf = and Pf = Maxwell's equations give rise to a differential K[A}-module M 
with generators ei, e 2 , e 3 , di, d 2 , d 3 , hi, h 2 , h 3 ,bi,b 2 , b 3 satisfying 

5 x di + 5 y d 2 + 5 z d 3 = = 5 x bi + 5 y b 2 + S z b 3 , 

5 y e 3 - 5 z e 2 + 5 t bi = = 5 y h 3 - 5 z h 2 - 8 t d x , 

S z ei - 5 x e 3 + 5 t b 2 = = 5 z hi - 5 x h 3 - 5 t d 2 , 

S x e 2 - d~ y ei + 5 t b 3 = = 5 x h 2 - 5 y hi - 8 t d 3 . 

Then M is isomorphic to the factor module of a free K[5 X , 5 y , 5 Z , 5t]-module with free gener- 
ators pi, . . . ,pi 2 by its submodule N generated by 

G = {S x p 7 + 5 y p 8 + S z p 9 , 5 x p w + S y pu + 5 z pi 2 , 5 y p 3 - 5 z p 2 + 5 t pio, 
5 y P6 - S z p 5 - S t p 7 , S z pi - 5 x p 3 + 5 t p u , 5 z p A - 5 x p 6 
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S X P2 - S y P! + S t p 12 , S x p 5 - S y p A - S t p 9 }. 

We define an admissible order -< by 

S a x *S^S a /S?e h -< «=► 
[a x + a y + a z + a t , ji, a x , a y , a z , a t ) 
<iex (b x + b y + b z + b t ,j 2 , b x , b y , b z , b t ). 

Then G is a Grobner basis and by Proposition IV. 6 the differential dimension polynomial 
associated with Maxwell's equations for vanishing free current density and free charge density 
is given by 

if s 1 4 19 3 55 2 137 

4 o 4 o 

Difference dimension polynomial for forward difference scheme 

Let K be an inversive difference field with basic set A = {8 X : x i— >■ x + 1, 8 y : y i— > 

?/ + 1, 5 Z : 2 i— >• z + 1, 5t : t h- >• t + 1}. If we replace every occurrence of Si in G by Si — 1, 

where i G {x, y,z, t}, we obtain a set 

G = {S x p 7 -p 7 + SyPs -ps + S z p 9 - p 9 , 



SxPl 


-Pio - 


t- SyP! - 


- Pn + S z p 2 


~Pl2, 


SyP3 


-P3- 


S Z P2 - 


P2 + S t Pl ~ 


Pio, 


$yP& 


-Pe- 


SzPh - 


P5 ~ StPl ~ 


P7, 


SzPl 


-Pi- 


S x Ps - 


P3 + StPl ~ 


Pn, 


SzP4 


-P4- 


SxPe - 


P6 ~ StP8 ~ 


Pa, 


SxP2 


-P2- 


SyPl ~ 


Pi + S t p 2 - 


P12, 


SxPb 


-Pb- 


SyPA ~ 


Pa ~ S t p 9 - 


Pq}- 



associated with Maxwell's equations for vanishing free current density and free charge den- 
sity. Applying the Buchberger algorithm we obtain the following 80-element Grobner basis 
G of the associated K[a x ,a y ,a z , a t , j3 x , j3 y , j3 z , (3 t ]-submodule of the free module with free 
generators pi, . . . , p\ 2 . 

G = {(3 x (3 y (3 z pi 2 + /3 x /3 y /3 z pu - (3 x (3 y p 12 + (3 x (3 y (3 z p w - /3 x f3 z p n - (3 y (3 z p w , (3 y (3 z pi 2 - a x (3 y (3 z p w 

+ PyPzPll ~ (3yPl2 + PyfizPlO ~ PzPll, PxPyPzPd + PxPyPzPS ~ PxPyPd + fixfiyfizPl 
- PxPzVi ~ PyPzPTi /3y/3zP9 ~ ax/3y/3 z P7 + PyPzPs - P>yP% + PyPzPl ~ PzP8, PyPzPtPl 
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- PyPzPtPG ~ PyPzP7 + PyPzPtPb + PzPtPd ~ PyPtPb: ~PyPtP7 + ^zPyPtPh + PyPtPd + PyP7 
~ PyPtPb ~ PtP6, ~PxPzPtPs + PxPzPs ~ PxPzPtPd + PzPtP% + PxPzPtPi ~ PxPtPl PxPtPs 
~ PxP8 + PxPtPQ + Oi z j3 x f3 t P4 ~ PtP6 - PxPtPl, PxPyPtP9 ~ PxPyP9 ~ PxPyPtPb + PxPyPtPA 

+ PyPtPs - PxPtPl PxPtP9 ~ PxP9 ~ Oy PxP&i ~ PxPtP5 + PxPtPl + PtP5, 'PyPzPtPlO 
+ PyPzPlO - PyPzPtPz + PyPzPtP2 + PzPtPl ~ PyPtP2, PyPtPlO ~ PyPlO + ^zPyPtPl + PyPtPz 
~ PyPtP2 ~ to, PxPzPtPll ~ PxPzPll ~ PxPzPtPz + PzPtPz + PxPzPtPl ~ PxPtPl, ~PxPtPu 
+ PxPll + PxPtP-A + CtzPxPtPl ~ PtPz - PxPtPl, -PxPyPtPia + PxPyP\2 ~ PxPyPtP2 
+ PxPyPtPl + PyPtP2 ~ PxPtPl, ~PxPtPl2 + PxP\2 ~ CtyPxPtPl ~ PxPtP2 + PxPtPl + PtP2, 

- PtP9 +P9~ a x p t p 5 + oi y p t j) A + p t p 5 - p t p 4 , -PtPs + a x PtPe +Ps~ PtPe ~ &zPtPi 
+ PtP4, -a y PtPs - PtPi + a z PtP$ + PtPe + Pi ~ PtPs, ~PtP\2 + P\2 + a x p t p2 - a y PtPi 

- PtP2 + PtPl, -PtPll + Pll - OixPtPZ + PtP3 + &zPtPl - PtPl, -PtPlO + PlO + UyPtPz 

- Oi z PtP2 - PtP3 + PtP2, OtyPzPn - PzP\2 + UxPzPlO - PzPll + Pl2 ~ PzPlO, CtyPzPfi ~ PzP9 
+ OLxPzPl ~ PzPS +P9- PzP7, -OCtP9 +P9 + OixPb ~ &yPi ~ £>5 + PA, -«iP8 + PS ~ &xP% 

+ p 6 + a z pi - p 4 , -a t p 7 + a y p 6 + p 7 - a z p 5 -p 6 +Ps, a t p u - P12 + a x p 2 - a y p x - p 2 
+ Pi, ct t p n - p u - a x p 3 +p 3 + a z pi - p 1 , a t p w - p w + a y p 3 - a z p 2 - p 3 + p 2 , a z p 12 
+ OiyPn - pyi + a x p 1Q - p n - p 10 , a z p 9 + a y p s - p 9 + a x p 7 - p 8 - P7} 
UiaiPiPj-pj I % G {x,y,z,t},j e {1, . . . , 12}} . 

Applying Proposition IV. 6 we obtain that the difference dimension polynomial associated 
with the forward difference scheme is given by 

0(r) = 4r 4 + 18r 3 + 35r 2 + 31r + 12. 

Difference dimension polynomial for symmetric difference scheme 

Let K be an inversive difference field with basic set A = {S x : x \— > x + l,8 y : y y y + 1,S Z : 
z \ — y z+1, 5t '■ 1 1— > t+1}. Using the symmetric difference scheme we replace every occurrence 
of Si in G by — S^ 1 ) (i G {x, y, z, t}) and arrive at a set 

G = {S x p 7 -P7 + S y Ps ~PS + $zP9 ~ P9, 

SxPio ~ Pw + SyPu - pn + S z p u - P12, 

SyP3 -P3- $zP2 -P2 + S t p 10 ~ p W , 
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SyPe -pe- $zP5 -Pb- $tP7 - P7, 



S z pi - Pi - S x p 3 -p 3 + S t pu - Pu, 



bzPi - Pa- S x p 6 - pe- $tPs - Ps, 



8xP2 - Pi- SyPi -Pi + 5 t p 12 - p 12 , 

SxP5 ~Pb- 8yP4 ~Pi- 5tP9 - Pq}- 



Proceeding as above we obtain that the corresponding difference dimension polynomial is 
given by 

56 

(j){r) = 4r 4 + — r 3 + 36r 2 + 4r + 22. 
3 

Comparing difference dimension polynomials computed for the forward and symmetric 
difference schemes we can conclude that the strength of the system of difference equations 
obtained via forward difference scheme is higher than the strength of the system obtained 
with the use of symmetric difference scheme. This time we obtain that the forward scheme 
is more preferable from the point of view of strength. 

Example IV. 9 (Equations for electromagnetic field given by potential) 

An electromagnetic field can be defined by the differential equations describing its potential, 
c/ 10 (Ex. 9.2.6.). The corresponding system of PDEs, which involves four unknown functions 
ijj\{x\, . . . , x 4 ), . . . , i>4(xi, . . . , X4), is as follows. 



Differential dimension polynomial 

d 

Let K be a differential field with basic set A = {5i = — — | i = 1, . . . ,4}. Then equations 
(7) and (6) give rise to a differential K[A]-module M with generators mi, . . . ,m 4 satisfying 
for i — 1, ... ,4 the defining equations 







1 



4 
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Then M is isomorphic to the factor module of a free K[A]-module with free generators 
ei, . . . , e4 by its submodule N generated by 



.3=1 



3=1 



(8) 



Defining an admissible order -< by 



(ai + a 2 + a 3 + 04, ji, «i, ^2, «3, «t) 

<lcx (&1 + &2 + &3 + 64, J2, 6l, &2 > fo 3, 64), 

we obtain the following Grobner basis G for N. 

G = {5 2 e 3 + 5 2 2 e 3 + 5le 3 + 5 2 3 e 3 
5fe 2 + 5 2 2 e 2 + 5\e 2 + 5\e 2 
5\e + 5\e + 5\e + 5\e 
<?ie + 5 2 e 2 + 5 3 e 3 + 5 4 e 4 
<5 2 e 4 - 5i<5 4 e + <5 2 e 4 - <5 2 <5 4 e 2 + <5 2 e 4 - <5 3 5 4 e 3 }. 

Applying Proposition IV. 6 we obtain that the differential dimension polynomial associated 
with (7) and (6) is given by 

3 11 2 17 
0(r) = r d + yr 2 + yr + 4. 

Difference dimension polynomial for forward difference scheme 

Let K be an inversive difference field with basic set A = {5i : x\ 1— > x\ + 1 | % — 1, . . . , 4}. 

Replacing every occurrence of 5k in (8) by 5k — 1 (k = 1, ... ,4) we obtain that the desired 

dimension polynomial is the A* -dimension polynomial associated with the factor module of 

4 

the free K[A*]-module E = K[A*]e{ by its K[A*]- submodule generated by the set 

i=i 



e 3 - e 3 > U 



.i=i 
4 

(<5|ej — 25jej + ti — 5i5jtj + <5j&j + SjSj — e^) 



i = l 4V. 
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Exploring the idea described at the end of subsection A, let us consider K as a difference 
field with basic set 

£ = {a,i : Xi H> Xi + 1, f3i : x { H> x { - 1 | i = 1, . . . , 4} 
and let the admissible order -< be given by 

:<==^ (a 1 + b 1 -\ h a 4 + 6 4 , z, ai, . . . , a 4 , 6i, . . . , 6 4 ) 

<icx (ci + d 1 H h c 4 + d 4 , i, Ci, . . . , C4, C?i, . . . , C?4j. 

Using the Maple package "Ore-Algebra" 21 for computing a Grobner basis of the K\£?[- 
submodule generated by 

X/ 1 ." / - j u 

4 

^^^^^ ^Qjj 6^ ^iOlj€i-t^ I Olq^OljGj I Glr^Qj I Oly Cj 6j^j 



i = 1. 



we obtain the set of leading terms of the Grobner basis 

{a 4 /3 8 ei, a 3 /3 7 ei, a 2 /3 6 ei, aiA^i, afei, a 4 /3 8 e 2 , ajj&ei, 
aiPefrPle!, a 3 (3 7 e 2 , a 2 (3 6 e 2 , ai/3 5 e 2 , a 2 e 2 , a 4 /3 8 e 3 , 
/3 8 e 4 , anfePefae!, al/3 5 f3 6 e u a 2 (3%(3 7 (3%e 1 , al(3 5 e 2 , 
aif3 6 (3 7 f3ie 2 , a 3 (3 7 e 3 , a 2 (3 6 e 3 , ai(3 5 e 3 , a 2 e 3 , a 4 e 4 , a 3 /3f/3f/3fei, 
a^(5 & (5 7 e 2 , al(3 5 (3 6 e 2 , a 2 (3l(3 7 (3le 2 , a%(3 5 e 3 , ai/3 6 /3 7 /3|e 3 , a 3 (3 7 e 4 , 
a 2 f3 6 e4, ai/3 5 e 4 , a 2 e 4 , /3|/3e/3f/3|ei, a 3 /3 2 /3|/3fe 2 , a^AsA)/^, 
tt 3 /35/3ee 3 , a 2 /3 5 2 /3 7 /3 8 2 e 3 , a 2 /3 5 e 4 , aiAj/^, Z 3 !/ 3 !/ 5 ?/ 3 !^, 

Applying Proposition IV. 6 we compute the difference dimension polynomial associated with 
a forward difference scheme for (7) and (6) to fre 

0(r) = 15r 3 - ^r 2 + yr 2 + 2. 

Difference dimension polynomial for symmetric difference scheme 
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Let K be an inversive difference field with basic set A = {5i : Xi > x,i + 1 | i — 1, . . . ,4}. 
Replacing every occurrence of 5k in (8) by \(5k ~ ^/T 1 ) (k = 1, ... ,4) we obtain the set 



SiSjCj + 5 i 5j 1 e j + S^SjCj - S^S^ej) 



which generates a K[A*]-submodule N of the free K[A*]-module M = i^[A*]e« such 

i=l 

that the difference dimension polynomial of our system of difference equations is the A*- 
dimension polynomial of M/N . 

Using the approach described at the end of subsection A, we treat K as a difference field 
with basic set 

£ = {ctj : Xi t-> Xi + 1, (3i : Xi H- Xj — 1 | i = 1, . . . , 4} 
and consider the admissible order -< given by 

a^aZafatfr&^^e, -< a? a? a? a? fl 1 ^^^e 3 
:^=> (ai + bi + • • • + a4 + 64, i, a±, . . . , 04, b±, . . . , 64) 

<icx (ci + di H h c 4 + d 4 , i, ci, . . . , c 4 , di, . . . , d A ). 

Once again using the Maple package "Ore_Algebra" for computing a Grobner basis of the 
K[Tj\-submodule generated by 



U 



f 4 1 



'>/'>/., + 0, , + h>\ j< , - PiPjej) 



we obtain the set of leading terms of the Grobner basis 

{a 4 /3 8 ei, a 3 /3 7 ei, a 2j 5 6 ei, ai/3 5 ei, afei, a 4 /3 8 e 2 , a|/3 5 ei, a 3 /3 7 e 2 , 
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a 2 (3 6 e 2l aci/3 5 e 2 , ocfe 2 , a 4 /3 8 e 3 , /3|e 4 , ap 5 (3 6 e u al/3 5 e 2 , a 3 (3 7 e 3l 
a 2 (3 6 e 3 , ai(3 5 e 3 , afe 3 , a 4 e 4 , /3%f3 6 (3 7 f3 8 ei, al/3 5 (3 6 e 2 , 
al/3 5 e 3 , a 3 (3 7 e 4 , a 2 (3 6 e 4 , «i/3 5 e 4 , a\e A , f3lf3 6 (3 7 f3 8 e 2 , alf3 5 f3 6 f3 7 e 2 , 
"iA;/^, al/3 5 e 4 , /3f/3 6 /3 7 /3 8 e 3 , a^AsA^, a|/3 5 /3 6 e 4 , /3|/? 6 /5 7 e 4 } 

Allying Proposition IV. 6 we compute the difference dimension polynomial associated with 
a symmetric difference scheme for (7) and (6) to be 

0(r) = 16r 3 - 8r 2 + 24r + 8. 

Comparing difference dimension polynomials computed for the forward and symmetric dif- 
ference schemes we see that in this case, as in the previous example, the forward scheme is 
more preferable from the point of view of strength. 

CONCLUSION 

We have developed a method for evaluation of the strength of systems of partial differen- 
tial and difference equations based on the computation of the corresponding differential and 
difference dimension polynomials. We have also presented algorithms for such computation 
that extend the Grobner basis technique to the cases of differential, difference, and inversive 
difference modules. Finally, we have determined the strength of some fundamental systems 
of PDEs of mathematical physics and the strength of the corresponding systems of partial 
difference equations obtained by the forward and symmetric difference schemes. 
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